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We consider Bayesian constraints on standard isocurvature baryon models from the slope and normalization 
of the anisotropy power spectrum detected by the COBE DMR experiment in their two year maps. In 
conjunction with either the amplitude of matter fluctuations <ts or its slope, all open models are ruled 
out at greater than 95% confidence, whereas cosmological constant A dominated models are constrained 
to be highly ionized. By including the COBE FIRAS 95% confidence upper limit on spectral distortion, 
we further reduce the available parameter space for A models by excluding these highly ionized models. 
These constraints define a single class of barely viable standard models which makes definite and testable 
predictions for degree scale anisotropics and large scale structure. 
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Obtaining life is a matter of timeliness. 
Losing it is a matter of compliance. 
Repose in timeliness and dwell in compliance 
Then sorrow and joy can never enter. 

-Chuang-tzu 

1. Introduction 

The original baryon isocurvature scenario for structure formation (Peebles 1987a, b) presents a simple 
and attractive alternative to the standard cold dark matter (CDM) cosmogony. It simultaneously satisfies 
observations which require a low density universe fio — 0.2 — 0.3 (e.g. Dekel et al. 1993), forms structure 
without the aid of hypothetical dark matter, and can alter light element nucleosynthesis sufficiently to make 
an fio = baryonic universe acceptable (Gnedin & Ostriker 1992). Moreover recent measurements of a 
large Hubble constant Hq = lOOTikms"! Mpc'^ h = 0.80 ±0.17 (Freedman et al. 1994) would be easier to 
accommodate in such a low density universe. 

Unfortunately, when normalized to the COBE DMR detection (Smoot et al. 1992), the open universe 
manifestations of this model appear to be inconsistent with several observations of cosmic microwave back- 
ground anisotropics at degree (Chiba, Sugiyama, & Suto 1994; Hu & Sugiyama 1994, hereafter HS94), and 
arcminute scales (Efstathiou & Bond 1987; Hu, Scott, & Silk 1994). The model generically suffers from 
excess small scale power. However, given the present uncertain status of CMB anisotropy detections at 
degree to arcminute scales (see e.g. Wilkinson 1994), it is perhaps premature to rule out models on these 
grounds. 

The excess small scale power required in this scenario also appears as a steep slope in the large angle 
anisotropy spectrum (Sugiyama & Silk 1994). This prediction conflicts with the flat spectrum measured by 
the COBE DMR experiment (Gorski et al. 1994; Bunn, Scott, & White 1995). In this Letter, we quantify 
this constraint on the standard isocurvature baryon model by employing the techniques of Bunn & Sugiyama 
(1995) to analyze the two year COBE DMR maps. Unlike previous treatments (Chiba, Sugiyama, & Suto; 
HS94), we also use the full information in the COBE sky maps to determine the normalization as opposed 
to merely the rms fluctuation at 10 degrees. This causes a 10% boost in the amplitude of fluctuations in 
open models. However we further extend prior treatments by considering flat low fio, cosmological constant 
A models whose predictions are somewhat more in accord with observations. The boost in amplitude can 
be up to 30% in these models. The corresponding change in the matter fluctuation amplitude <ts is relevant 
for simulations of large scale structure formation. Finally, employing spectral distortion constraints from 
the COBE FIRAS experiment (Mather et al. 1994), we nearly close the parameter space available to these 
baryon isocurvature models. For the small class of models remaining, we present the predictions for degree 
scale anisotropics and large scale structure. 
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2. General Features 

In the standard baryon isocurvature model, the universe consists of photons, baryons, and three families 
of massless neutrinos only. Initial entropy perturbation, i.e. fluctuations in the baryon-photon and baryon- 
neutrino number densities, are assumed to take the form of a pure power law in k, cx k" where the 

wavenumber k is related to the eigenvalue of the Laplacian k as k^ = k^ + K, with K = —Hq(1 — — ^^a) as 
the curvature (Wilson 1983). Here is the fraction of the critical density contributed by the cosmological 
constant. The A dominated models which we consider here are flat for simplicity, i.e. K = 0. In this case, 
k = k and represents an ordinary Fourier mode of the perturbation. This limit is also appropriate for large 
scale structure measurements which are affected by perturbations on scales that are well under the curvature 
radius. 

Since there is no ah initio mechanism for generating the entropy perturbations, the index n is fixed 
by measurements of large scale structure today. Isocurvature perturbations evolve such that below the 
photon diffusion scale, the initial entropy fluctuations become the density perturbations that seed large scale 
structure. This implies that the observational constraints of an P{k) 2^ k~^ power spectrum at large scale 
structure scales (e.g. Peacock & Dodds 1994) implies an n ~ —1 initial power law in the model. Numerical 
simulations which take into account non-linear modifications of this picture confirm this result (Suginohara 
& Suto 1992). At the largest scales however, isocurvature conditions prevent the formation of density 
perturbations leading to a steep P(k) cx (k'^ — AK)'^k" , i.e. an n + 4 power spectrum below the curvature 
scale. This sharply rising spectrum of fluctuations on scales relevant for CMB anisotropics is particularly 
dangerous when normalized at the large scales by the COBE DMR measurement. 

It may thus seem that the model can be ruled out by merely considering the implied amplitude of the 
matter power spectrum at the 8/}~^Mpc scale. In an unbiased scenario of galaxy formation, which is expected 
in these baryon only models (Cen, Ostriker, & Peebles 1994), observations require <ts — 1. However, the 
baryon isocurvature model has an additional degree of freedom to save it. Since Silk damping (Silk 1968) 
does not destroy entropy fluctuations, the large amount of small scale power in the model allows for collapse 
of objects immediately following recombination. This could lead to sufficient energy input to reionize the 
universe (Peebles 1987). Because Compton drag prevents the growth of structure, the ionization history 
can be tuned to provide the right ratio of matter to temperature fluctuations. Following Gnedin & Ostriker 
(1992), we assume that a fraction Xg of the electrons were reionized at z ~ 800. For complications due to a 
multi-staged ionization history and compact baryonic object formation, see HS94. 

Reionization also leads to significant and observable consequences for the CMB. Large primary fluctu- 
ations from the acoustic oscillation phase (see e.g. Hu & Sugiyama 1995) are exponentially damped with 
optical depth below the horizon at the new last scattering surface. Secondary anisotropics are generated 
due to Doppler shifts off of moving electrons at last scattering. These are damped under the thickness of 
the last scattering surface due to cancellation of redshifts and blueshifts as the photon travels across many 
wavelengths of the perturbation. Thus it is almost always the case that higher ionization implies smaller 
anisotropics under the angle that the horizon subtends at last scattering. We plot the anisotropics in a A 
model as a function oix^ in Fig. 1, where the rms anisotropy is related to via (| AT/Tp) = '^{2l-\-l)CilAi:, 
with I as the multipole number of the spherical harmonic decomposition of anisotropics on the sky. Open 
universe examples are displayed in HS94. 
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Figure 1. COBE normalized anisotropics in the A model as a function of ionization. As 
the ionization level increases from = to 1 as listed in the figure, the damping reaches 
to larger angles making the COBE slope shallower. Open universe models suffer less from 
this effect at large angles due to geodesic deviation. Fluctuations are also regenerated on 
the new last scattering surface. With high enough ionization, this can once again steepen 
the COBE slope. The COBE normalization also sets the level of matter fluctuations <ts for 
a fixed thermal history. The ratio r of amplitudes between the more complete likelihood 
analysis used here and the 10° rms normalization is shown. The most promising model which 
currently escapes constraints from the COBE slope, normalization, and spectral distortion 
measurements is shown (thick line, Xg = 0.1). This should be compared with the standard 
CDM {h = 0.5, fift = 0.05, Qc = 19.9yuK) model (dashed line). Measurements in the range 
£ ~ 20 — 200 can further help to distinguish the models. 

One exception to this damping rule is the second order Doppler contributions from the Vishniac effect 
(Ostriker & Vishniac 1986; Vishniac 1987) which is not included in Fig. 1. This effect is uncovered at 
arcminute scales, where other first order effects have suffered severe thickness damping, and is extremely 
sensitive to the amplitude of the matter perturbations. It thus is only important for highly ionized, late last 
scattering scenarios (Hu, Scott, & Silk 1994; HS94). 

Finally ionization also implies that the electrons have been heated to a temperature above that of the 
CMB. This implies that Compton scattering will lead to spectral distortions in the CMB as photons are 
upscattered in frequency by the electrons. The distortion is described by the Compton-j/ parameter defined 
as y = J drkiTe — T)jmf.c'' where Tg and T are the electron and CMB temperatures respectively and r is the 
optical depth to Compton scattering. Explicitly dr = (xene(JTc)dt, where Xe is the ionization fraction, rig is 
the electron number density, and ax is the Thomson cross section. Thus the greater the level of ionization, 
the larger the spectral distortion in the CMB. 
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3. Model Constraints 

When extended to large scales, the steep initial spectrum required by large scale structure conflicts 
with the flat anisotropy spectrum measured by the COBE DMR experiment. Without reionization, the 
predicted COBE DMR slope is approximately rieff — 2 (Sugiyama & Silk 1994) compared with observational 
constraints of rieff = 1-31q'37 (with quadrupole, Bunn, Scott, & White 1995). The fact that rieff is only 
weakly dependent on n and is somewhat shallower than one might expect from the n + 4 behavior of the 
matter power spectrum is discussed in Hu & Sugiyama (1995). 

Reionization tends to suppress small angle anisotropics and can mitigate a steep initial spectrum. How- 
ever, if the ionization is too great, secondary anisotropics generated on the last scattering surface will counter 
this effect (see Fig. 1). Furthermore, reionization has no effect for angles much larger than that subtended by 
the horizon at last scattering. At the COBE scale, open models will thus be less affected by reionization than 
A models, since geodesic deviation carries the same physical scale at last scattering to a much smaller angle 
on the sky today. Lesser effects can be attributed to raising the baryon content through Q.},h^ which delays 
last scattering and increases the physical scale of the horizon. However even for flat models, projection effects 
due to the distance to the last scattering surface depend strongly on fio and counters the fi;, dependence 
in these fio = baryonic models. Furthermore, the late integrated Sachs- Wolfe effect (Sachs- Wolfe 1967; 
Hu & Sugiyama 1995) boosts the low order multipoles slightly as fio decreases. In the range of interest, 
decreasing fio leads to a shallower COBE slope. High Xg, high h, low fio, A models therefore offer the best 
prospects of bringing down the COBE slope. 

Employing the two year COBE DMR sky maps and the analysis methods developed by one of us (EB), 
we can place an upper limit on the primordial index n for open and A models fixed by fio, h, and Xg. 
We expand the two-year DMR data in a set of basis functions which are optimized to have the maximum 
rejection power for incorrect models (Bunn & Sugiyama 1995; Bunn, Scott, & White 1995). We then use 
the 400 most significant terms in this expansion to compute the likelihood functions for a variety of models. 
To set limits on n and the normalization Q, the rms quadrupole, we assume a prior distribution which is 
uniform for all Q and n < 0. Spectra with n > are unphysical due to non-linear effects which regenerate an 
n = 0, P(fc) cx fc'* large scale tail to the fluctuations (Zel'dovich 1965; Peebles 1980). The constraint in the 
crucial n ~ — 1 regime is not sensitive to the details of this cutoff. Shown in Fig. 2 are the 95% confidence 
upper limits imposed on n by integrating over the normalization Q to form the marginal likelihood in n. As 
expected, all open models with n ~ — 1 are ruled out regardless of ionization fraction, whereas highly ionized 
A models remain acceptable. 

With the maximum likelihood value for the normalization Qc of the model, we predict the amplitude 
of matter fluctuations as- The likelihood value for the normalization tends to boost the amplitude over the 
COBE DMR 10° rms normalization value of 30//K (Bennett et al. 1994) by a factor r = QclQw — 1-1 for 
open models and low ionization A models. The difference is more significant in highly ionized A models due 
to the damping of the anisotropy spectrum. The boost is on the order r ~ 1.3 for a fully ionized A model 
(see Fig. 1). This effect appears also in the CDM model with a greater magnitude in fact. The effect of 
the low quadrupole in the data on the 10° measure (Bunn, Scott, & White 1995) artificially supresses the 
amplitude. In all cases, the likelihood analysis provides the better normalization by including the full data 
set and minimizing the effects of cosmic variance. In Fig. 2, we thus plot the value of ag corresponding 
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Figure 2. Constraints on the primordial spectral index n and ionization fraction x^. The 
COBE DMR slope imposes a 95% upper confidence limit on n which is weakened as the 
ionization fraction increases, due to damping of the primary fluctuations, until a turning 
point at which fluctuations are significantly regenerated by the Doppler effect on the new 
last scattering surface. The COBE FIRAS constraint on spectral distortions through the 
Compton-j/ parameter sets an upper limit on the ionization fraction. Here a conservative 
Te = 5000K is assumed. The more realistic Tg = lOOOOK is also shown (dotted lines). 
The COBE DMR normalization also sets the level of matter fluctuations at the 8 h~^Mpc 
scale (78. (a) No open model simultaneously satisfies all the observational constraints, (b) 
For A models, a small region of parameter space is open for high h, low fio models. The 
full anisotropy spectrum for the most promising model fio = 0.2, h = 0.8, n = —1.15 and 
Xe = 0.1 is displayed in Fig. 1 and the matter power spectrum in Fig. 3. Even this model is 
ruled out with the more realistic Tg. 

to this normalization as a function of ionization history and spectral index. The suppression of fluctuation 
growth in a highly ionized universe must be compensated by a steeper spectral index n. Notice that even 
ignoring limits on the large scale structure slope, all open models which satisfy the COBE slope are ruled 
out. 

Even though highly ionized A models can survive constraints on the COBE slope and the large scale 
structure normalization, they run into difficulties with the low upper limit on spectral distortions imposed by 
the COBE FIRAS experiment, y < 2.5 x 10~® (95% CL). If the intergalactic medium is collisionally ionized, 
the electron temperature must be Te > 10, OOOK (see e.g. Gnedin & Ostriker 1992). To be conservative, we 
take Te = 5000K. The corresponding limit from the Compton-j/ parameter may be avoided by more exotic 
ionization schemes which attempt to inject as little energy as possible into the electrons (e.g. neutrinos 
decaying via a 13.6 eV photon). However since we generically expect at least a few eV excess energy above 
the ionization threshold, 5000K is a reasonable minimal estimate of the electron temperature. Calculations 
by Tegmark & Silk (1994) of photo-ionized models, which include feedback from the CMB through Compton 
cooling, support this conclusion. With this constraint, even A models fall from favor. Only high h models 
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Figure 3. Power spectrum k^P(k) for an allowed A model (fio = 0.2, h = 0.8, Xe = 0.1, 
n = —1.15). For comparison, a parameterized F = 0.2 CDM model (Efstathiou, Bond, & 
White 1992) is shown normalized to (Xg = 1 which is known to fit the shape of the large 
scale structure data at 10"'^ < k/h < lMpc~^. To facilitate more detailed comparisons, 
we have also provided a simple fitting formula involving the isocurvature transfer function 
P(k) oc |T(fc)5(fc)p and the maximalJeans scale kj which is motivated by the perturbation 
analysis of Hu & Sugiyama (1995). 

have a small window of parameter space open, which in fact closes if a more realistic Tg = lOOOOK is assumed. 
For these early ionized scenarios, this constraint largely obviates the need to impose limits from the Vishniac 
effect (HS94). 

Finally, it is worthwhile to mention that once the observational situation at degree scales settles down, 
one can at the very least impose a lower limit on the ionization fraction (see Fig. 1 and HS94). Significant 
reionization (xe ^ 0.01) is necessary in these models to avoid large degree and arcminute scale fluctuations in 
the CMB. Moreover, with the rapidly increasing number of experiments, the sample variance associated with 
the measurements (Scott, Srednicki, & White 1994) will decrease to the point where all baryon isocurvature 
models can be distinguished from the CDM model. In Fig. 1, we have plotted the predictions for the most 
promising baryon isocurvature model (thick solid line) in comparison to the CDM model (dashed line). The 
shape to the rise of the prominent peak around £ ~ 100 — 200 and the fall off thereafter may be used to 
help distinguish the models. In Fig. 3, we plot the matter power spectrum for the same model and include 
a simple fitting formula which may facilitate the comparison with large scale structure measurements. As 
a simple comparison, we also show a F = 0.2 CDM model which is known to fit the slope of the observed 
spectrum well (Efstathiou, Bond, & White 1992). These models may consequently run into problems with 
the large and small scale regimes of large scale structure measurements as well as indications of a smooth 
power spectrum (Peacock & Dodds 1994). With these additional constraints, one may hope to close off the 
already small window of parameter space available to the model. 
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4. Discussion 

standard baryon isocurvature models generically run into conflict with CMB observations even at COBE 
scales. Open models are ruled out by a combination of the COBE spectral slope and implications of the 
normalization for the matter power spectrum. Whereas these two considerations leave a large window of 
acceptable A models, the inclusion of the COBE FIRAS constraint on spectral distortions even in a relatively 
conservative fashion is sufficient to drastically reduce the available parameter space such that a tuning of 
the ionization history, fio and h must be involved. 

With the current state of affairs in which none of the simplest models for structure formation fare well 
in comparison with all the observations of the CMB and large scale structure, it is perhaps unwise to dismiss 
this scenario as entirely unviable. The general idea of isocurvature seeded fluctuations may of course be 
saved by introducing more free parameters. 

The original model employs two simplifying assumptions: a power law initial spectrum and a constant 
ionization fraction after reionization. Since open models run into difficulties by predicting a steep COBE 
slope, the former assumption must be dropped to save them. In fact, for the open models there is some reason 
to believe that the spectrum may possess non-trivial structure at the curvature scale (Lyth & Stewart 1990; 
Ratra & Peebles 1994; Bucher, Goldhaber, & Turok 1994). Note however that unlike the open adiabatic 
case, power law behavior in gravitational potential fluctuations is equivalent to power law behavior in the 
entropy fluctuation (Hu & Sugiyama 1995), which serves to eliminate a potential ambiguity of the open 
model. Moreover, tentative indications of a break to a steep rise in anisotropics at degree scales (Scott & 
White 1994) may require further adjustment of the initial spectrum even below the curvature scale in the 
open model. 

On the other hand, more complicated ionization histories and compact baryonic object formation can 
be employed to help design a more favorable A, but not open, model. The ionization history can fixed such 
that the relative normalization of the matter and radiation yields (jg = 1 (see HS94). However, since even 
the maximally damped open models violate the COBE slope constraint if n ~ —1, no tuning of thermal 
histories alone can save the open baryon isocurvature scenario. For A models, thermal history effects can 
also be employed to escape the COBE FIRAS constraints without giving up the damping benefits of a highly 
ionized model. This is because spectral distortions are a function of the total optical depth, whereas the 
damping of anisotropics is determined at last scattering where the optical depth equals unity. Thus late 
ionized scenarios may be more favorable. Unfortunately, the large amount of small scale power may make 
delayed reionization impossible. 

More radical solutions have also been proposed. Peebles (1994) suggests the addition of cold dark matter 
or defects, which may also provide sufficient freedom to save the model. Small admixtures of adiabatic 
fiuctuations may also be added. Yet with these ad hoc patches on the model, the appeal of the baryon 
isocurvature scenario is greatly reduced. 
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